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We derive the Josephson relation for a dilute Bose gas in the framework of an auxiliary-held 
resummation of the theory in terms of the normal- and anomalous-density condensates. The mean- 
held phase diagram of this theory features two critical temperatures, T^ and T* , associated with 
the presence in the system of the Bose-Einstein condensate and superfluid state, respectively. In 
this context, the Josephson relation shows that the superfluid density is related to a second order 
parameter, the square of the anomalous-density condensate. This is in contrast with the corre- 
sponding result in the Bose gas theory without an anomalous condensate, which predicts that the 
superfluid density is proportional to the BEC condensate density. Our findings are consistent with 
the prediction that in the temperature range between Tc and T* a fraction of the system is in the 
superfluid state in the absence of the BEC condensate. This situation is similar to the case of dilute 
Fermi gases, where the superfluid density is proportional to the square of the gap parameter. The 
Josephson relation relies on the existence of zero energy and momentum excitations showing the 
intimate relationship between superfluidity and the Goldstone theorem. 
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I. INTRODUCTION 

Superfluidity represents a fundamental property of 
matter that causes liquids to lose all their resistance to 
internal movement [1]. In the limit of dilute Bose gases, 
Josephson showed that the superfluid density[2] is pro- 
portional to the condensate density given by the square 
of the broken-symmetry order parameter, (j). The Joseph- 
son relation reads [2, 3] 
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lim - 
9">0 g2Gii(q,0) 



(1.1) 



where Gn (fc, 0) is the Fourier transform of the connected 
Green function corresponding to the expectation value 
( <j)* {x)(j){x') ), evaluated at zero energy transfer and mo- 
mentum transfer, q. The Josephson relation was re- 
derived in the framework of diagrammatic perturbation 
theory by Holzmann and Baym [4]. In both approaches, 
the Josephson relation implies that the Bose-Einstein 
condensate (BEC) and the superfluid state are intimately 
linked. Hence, the critical temperature, correspond- 
ing to the temperature where the BEC condensate flrst 
appears in the system, coincides also with the temper- 
ature where superfluidity sets in. Unfortunately, it is 
well-known that dilute Bose gases theories involving a 
single order parameter predict a flrst-order phase tran- 
sition, contrary to the second-order transition expected 
for the universality class of the Bose gas [5] . 
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Recently, we reformulated the theory of dilute Bose 
gases by treating the normal and anomalous densities 
on equal footing [6]. In the leading-order approximation, 
our auxiliary-fleld loop expansion predicts a second-order 
BEC phase transition. The leading-order auxiliary-fleld 
(LOAF) approximation is a conserving and gapless ap- 
proximation with a non-perturbative character and can 
be applied outside the regime of weakly-interacting parti- 
cles, while still preserving the correct weak-coupling lim- 
its [7]. LOAF produces the same slope of the linear de- 
parture of the critical temperature from the noninteract- 
ing limit derived by Baym et al. [8] using a large-N expan- 
sion for the critical theory. Unlike the large-N expansions 
developed by Baym et al. [8], the LOAF approximation 
can be used at all temperatures. 

Unlike all other theories of interacting dilute Bose 
gases, however, the LOAF phase diagram features two 
critical temperatures: the critical temperature T^, where 
the BEC condensate appears flrst, and a tempera- 
ture T* > Tc, where the anomalous-density conden- 
sate, A, turns on. We used the Landau two-fluid phe- 
nomenological model [9] to show [10] that in LOAF a 
superfluid state is present below T* and that the super- 
fluid density is proportional to [10]. Moreover, LOAF 
indicates that in the temperature range between and 
T* the system supports zero energy and momentum ex- 
citations, which correspond to composite-fleld Goldstone 
states. This is in addition to the usual Goldstone theorem 
for T < Tc, where the natural U(l) symmetry breaking 
of the theory leads to the usual BEC condensate as the 
Goldstone state below Tc. Introducing a hypothetical 
U(l) gauge vector meson into the system, we were able 
to show that the mass generated through the Anderson- 
Higgs mechanism ]11] can be related to the superfluid 
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density via the Meissner effect, which also gives that ps 
is proportional to . Hence, LOAF predicts that part 
of the system is in a superfluid state, in the absence of a 
BEG condensate, in the temperature range between the 
BEG critical temperature and the temperature where su- 
perfluidity sets in. According to Eq. (1.1), the LOAF ap- 
proximation appears to contradict the Josephson relation 
for an interacting dilute Bose gas. 

In this paper, we will derive the form of the Josephson 
relation in the LOAF model. We will show that in our 
auxiliary-field formalism, the superfluid density is indeed 
proportional to the square of the anomalous-density con- 
densate. The latter is identifled as a second order param- 
eter in LOAF. The Josephson relation for the superfluid 
density in LOAF reads 



8m2 
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(1.2) 



Here, ^'^^i is the Fourier transform of the connected prop- 
agator corresponding to the expectation value {A* A). 
Furthermore, we will study in detail the connection be- 
tween the Josephson relation and the Goldstone theorem 
in the LOAF model. 

The form of the Josephson relation described in 
Eq. (1.2) is not surprising. In fact, the Josephson re- 
lation for the case of interacting dilute Fermi gases was 
derived by Taylor and has the form [12] 
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where Dn is the Fourier transform of the connected prop- 
agator corresponding to the expectation value ( A*A ) 
for quasi-particle bosons of mass, uib, constructed as 
pairs of mass mp fermions (m^ = 2mp). Given that 
the gap parameter. A, is in effect an auxiliary fleld, 
A = (ip^Tp^), the propagator Dn is an auxiliary-fleld 
propagator, similarly to what is described in Eq, (1.2). 
It is important to note here that the fermionic counter- 
part of the LOAF model for dilute Bose gases we dis- 
cuss in this paper [13], is precisely the standard Bardeen- 
Gooper-Schrieffer (BGS) ansatz [14, 15] that underlines 
the derivation of Eq. (1.3). The fact that in LOAF the 
Josephson relation (1.2) is different from the classical ex- 
pression (1.1) is due to the fact that in our auxiliary fleld 
formalism, we treat the normal and anomalous density 
condensate condensates on equal footing. In this sense, 
the auxiliary-fleld formalism discussed here provides a 
unifled approach to the study of fermionic and bosonic 
atom gases. The detailed derivation of the LOAF ap- 
proximation for the case of dilute Bose and Fermi atomic 
gases can be found elsewhere ]13, 16]. 

This paper is organized as follows: We begin by re- 
viewing the calculation of thermodynamic variables in 
the auxiliary-fleld formalism in Sec. II. In Sec. HI we 
discuss the Goldstone theorem in the auxiliary-fleld for- 
malism, and we show that Goldstone modes are to be 



expected both for temperatures below Tc, but also in the 
temperature range Tc < T < T* . Turning our attention 
to the study of the superfluid phase, we begin by deriv- 
ing the superfluid density in the Landau two-fluid model, 
and we show explicitly that in LOAF the superfluid den- 
sity is proportional the square of the anomalous-density 
condensate in Sec. IV. In Sec. V, we calculate superfluid 
density in our microscopic theory and show the connec- 
tion with the current-current correlation function. Here 
we follow an approach similar to that discussed by Tay- 
lor et al. in the context of interacting Fermi gases ]17]. 
Finally, in Sec. VI we derive the Josephson relation for 
superfluidity in our auxiliary-fleld formalism and show 
that the square of the anomalous-density condensate is 
the second order parameter in LOAF. 



II. THERMODYNAMICS OF THE BOSE 
SYSTEM IN THE AUXILIARY FIELD 
FORMALISM 

We use here a grand canonical ensemble to describe the 
system in equilibrium. In this case, the normalized ther- 
modynamic probability density matrix p for the grand 
canonical ensemble is obtained by minimizing the en- 
tropy, deflned by S/k^ = — Tr[pln[/5] ], such that the 
average energy E = Ti[ pH] and average number of par- 
ticles N = Tr[ pN] are held constant. The solution for 
the minimization is given by 
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-|3(H-^iN) 



(2.1) 



Z[T, p, V] = e-^ "[^^A^'^l ^ Tr[ e"^ i"-^^)' 



where we have set f3 — l/(fcBT) and ft[T,p,V] is the 
grand potential. In our auxiliary fleld formalism, we 
will compute the thermodynamic grand potential using 
a path integral formalism. Once we have the potential, 
we can flnd all the thermodynamical properties of the 
system. For instance, using the second law of thermo- 
dynamics, TdS[E,N,V] ^ dE ~ pdN + pdV, we flnd 



(2.2a) 
(2.2b) 



From the second law, the energy is given by 

E = n + TS + pN . (2.3) 

In the following, we review the auxiliary fleld formalism 
for dilute Bose gasses [6, 16] using the imaginary time 
formalism. 
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A. Path integrals 

The partition function is given by a path integral [18], 
Z[V,^l,P] JJ D0D0*e-^l'^''^''^'^'''l , (2.4) 
where the EucHdian action S[(f>, (f)*; V, ii, /3] is given by 

s[^,^*;V,^l,p]^ /[dx]/:[0,r;A^], (2.5) 



(2.6) 



and where we have introduced the notation, 



The EucHdian Lagrangian density is given by 



CI 



(2.7) 



= - [ 0* (x) dr 0(x) - ^(x) dr r (x) ]+n[(b,r■,^i], 

where the Haniiltonian density is given by 



2m 



ix)+^\<P{x)\\ (2., 



Auxihary fields x(a:) and A(x) are introduced by means 
of the Hubbard-Stratonovitch transformation Lagrangian 
density, which takes the form 

raux A ^1 --^ [x{x)-xV2\^{x)\^f 

+ ^\Aix)~Xq^Hx)\\ (2.9) 

which we add to Eq. (2.7), giving a Euclidian action of 
the form. 



S[J,<^>;V,fi. 
_ 1 
^ 2 



(2.10) 



[dx] [dx'] M^) G-^\{x,x')(f>\x) 
\A{x)\^-xH^) 



■dx] 
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Ja(x) <P"{x) 



where 



G-'\ix,x')^S{x,x') ^^If 



with 



=h + dr, /l^"^ =h-dr 



(2.11) 



(2.12) 
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+ V^xix) - fi 



and 



J„(x) $"(a;) = r (a;)0(x) + J(a;)0*(a;) 

+ s{x) x{x) + S* (x) A{x) + S{x) Aix) 



(2.13) 



Here we have added currents and introduced a two- 
component notation using Roman indices a,b,c, ■ ■ ■ for 
the fields (/)(a;) and <j)*{x) and currents j{x) and j*{x), 

r(x) = (0(a;),0*(x)), (2.14a) 

(l>a{x) = {(f)*{x),(j){x)) , 

j-(a;) = (j(a;),r(a:)), (2.14b) 

ja(2^) = {i*{x),j{x)) , 

for a — 1,2, and a three-component notation using Ro- 
man indices i,j,k,--- for the fields xi^)- ^(2^)1 s-iid 
A*{xl 

x\x) = {x{x),A{x),A*{x)) , (2.15a) 
X^{x) = {x{x),A*{x),A{x)) , 

S-^a;) = (s(a;),S'(a;),S'*(x)) , (2.15b) 
Si{x) = {s{x),S*{x),S{x)) , 

for i — 0,1,2. For convenience, we also define five- 
component fields with Greek indices <i>" (x) and currents 
J"(x), 

= ( r{x),x'{x) ) = ( 0, 0* , X, ^, ^* ) , (2.16a) 
$a = {Mx),X^{x)) = {<l>*,(t),x,A*,A) , 
J"(a;) - (r(x),5X^)) = (j,J*,s,^,^*) , (2.16b) 
-'"0(2;) = {ja{x),Si{x)) = s, S'*,5) , 

The metric which raises and lowers indices is given by 



/G 1 0\ 
'10000' 
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Vo 1 0/ 



(2.17) 



The partition function is now a functional of the currents 
and given by a path integral over all fields, which we write 
symbolically as 



Z[J;y,^,/3] =6-^'^'^'''^^' (2.18) 

with S[^,J;V,fi,p] now given by Eq. (2.10). The ther- 
modynamic average value of the fields are given by 



;$"(x); 



1 SZ[J] 



Z SJa{x) 



6W[J] 



j=o 



5Ja{x) 



J=0 



(2.19) 



evaluated at zero currents, and the connected two-point 
functions are given by the Green function 



Gafiix^x') 



S^W[J] 



6Ja{x) (5 J^(x') 



J=0 



(2.20) 
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B. Effective action 



evaluated at the stationary points. Here 



After introducing the auxihary fields, the action is now 
quadratic in the cjy&elds and can be integrated out. The 
partition function is then given by 

Z[J] = e-^'['^l =mJi)x e-^"«[X'^l/' , (2.21) 



where W[J] = /3f2[J] and where 



'S'off[x, J] 
1 

^ ^2 



(2.22) 



[da;; 



[dx] [dx']ja{x)G\[x]{x,x')j\x') 

\A{x)\^^xH^) 



2A 



S,ix)xHx) 



-Tt[HG-'[x]{x,x)Y_ 



The dimensionless parameter e in Eq. (2.21) allows us 
to count loops for the auxiliary-field propagators in the 
effective action. 

Next we expand the effective about the stationary 
points Xoix), defined by dScs[x, J]/^Xi{x) = 0, i.e 

^ = V2{\Mx)\^ + nG[xo]{x,x)]/2} ~ s{x) 



Aoix) 



^l{x) + G\[xo]ix,x) + 2S{x) 



(2.23) 



where '^□(x) is given by 



ro[xo,J]{x) = V^oixoK^) + J [dx']G[xo]\{x,x')f{x') , 

(2.24) 

where (y9g[xo](2;) is a solution of the homogenous equa- 
tion, 



G-'\[Xo]{x)v^'o[Xo]{x) = Q. 



(2.25) 



Eqs. (2.23) are called the "gap" equations. The fields 
Xo[J](a;) at the stationary points are functionals of all 
the currents J" (x) . Expanding the effective action about 
the stationary point, we find 



SMX,J]^ SMxo,J] (2.26) 
+ 1 JJ[dx] [dx']D;^'[xo]ix,x') 

xix\x)-xU^)){x'ix')-xiix')) + --- 



where D^j^[xo]ix,x') is given by the second-order deriva- 
tives 



n-ir u '^^ 'S'cff [x°] 



(2.27) 
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fii.j5{x,x') +nij[xo]{x,x') 
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r7„- = I 1/2 I , 
1/2 



(2.28) 



and nij[xo](2;, a;') is the polarization, given symbolically 

by 

n^^'ko] = l{~Mxo] o V'^xo] o Mxo] (2.29) 
+ Tr[G[xo]o^'oG[xo]o^^]}, 



with 



SXt 



(2.30) 



V^^^'Ixo] - V' o G[xo] o + o G[xo] o V' 



We perform the remaining gaussian path integral over 
the fields Xi, obtaining the result for the grand potential 



eW[J]=So + SMxa,J] 



(2.31) 



■da;]Tr[ln[Z?-i[Xo, J](a;,a;) 



where 5*0 is a normalization constant. The fields are given 
by the expansion. 



e$"[J](x) = -e 



SW[J] 



5Ja{x) 

J] (a;) + e J] (a;) 



(2.32) 



We calculate the order e corrections to the fields from 
(2.32), evaluated at zero currents. We work these out 
explicitly for each field to zeroth order and find 



0(o)'^(a;) = ^g[Xo](x) 

x(°)(a;)=Xo(:i;), 
A(o)(a;)=Ao(a;), 



(2.33a) 
(2.33b) 
(2.33c) 



where (/3q[xo](^) is given by Eq. (2.25) and Xo{x) and 
Aq{x) are given by Eqs. (2.23), evaluated at s{x) — 
S{x) = 0. Diagrams for the first order fields <i>(^^ "[J](a;) 
are given in Ref. 16. 

The grand potential r[<i>] as a functional of the fields <i> 
(rather than the currents J) is constructed by a Legendre 
transformation. 



er[$]=e j[dx]Ja{x)^°'{x) + eW[J] (2.34) 
= ^^ + lJJ [d^'] Mx) G-^ixTbix, x') /(x') 
+ /|dx]{^M^^x!(^+lTr[ln[G-Mx](.x,x)]] 



[\n[D-'[xo,J]ix,x)]]} 
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Here Tq is a normalization constant. Then the currents 
are given by 



,5r[$] 



(5<I>"(a;) 

and the inverse Green functions by 

(52r[$] 



(2.35) 



(2.36) 



C. Leading order auxiliary field (LOAF) 
approximation 

For uniform r-independent systems, <i>"(x) = are 
all constants. Let us first define the effective potential 
14ff[$] as the grand potential per unit four- volume as a 
functional of the fields. Then from (2.34) for uniform 
systems, 



vM^]^e^ = Vo + x'\^\'-l 



A6*^ + A* 6^] 



+ ^-^^^^ + ^Tr[ln[G-Mx](x,x)]] (2.37) 

Here we have set x' = %/2x~M- Expansion of the inverse 
Green function in a Fourier series gives 



1 



Tr[ln[G-Mx](a^,^)]] 



(2.38) 



;^/(g.EMdet[G-(k,n)]] 



2/3 J {2tt) 
d^k 1 



(27r)3 2/3 

312 /3 



(27r)3 
where 



ln[ 1 - 



2 7,2 



2m 



Inserting this result into (2.37) gives 



2A 4A ■ ./ (27r)3 L 2 ' /3 



(2.39) 



(2.40) 



(x'+A*)^ 



However this expression is not finite. Expanding Wfe 
about fc — >■ oo, we find 



Wfc = Cfe + x' 



\Al 
2ek 



(2.41) 



These three terms are responsible for the divergences in 
the integral in Eq. (2.40). Our method of regularization 
is to subtract these three terms from the integrand and 



replace the coupling constant, chemical potential, and 
the normalization constant by renormalized values, which 
gives 



4A« 



2 

d3fc 

(2^ 

}■ 



1 r 

2 



2 + A* 

i^fc - e/c - X + 



2Xr 
\Al 
2ek 



(2.42) 



In dilute atomic gasses, the renormalized coupling con- 
stant is related to the s-wave scattering length, ag, by 
Xr = Anh? aQ/m. The regularization method is discussed 
in Ref. 16. 

Evaluating yeff[*i'] at the minimum where 



. 



yields the three equations. 



X 

-A* 



^A 



X + f^B. 

2Xr 



d3fc 

(2^) 



^ r 



ik+x' 
2ujk 



[2n{l3uJk) + I] 



A 
X^ 



-A 



r d^k r 1 
y (2^1 20^;^ 



2n(/3cjfc) 



1 - 



(2.43) 

(2.44a) 
(2.44b) 

\}- 

(2.44c) 
2^}' 



where n{x) — l/(e^ — 1) is the free particle Bose number 
distribution function. The particle density p is given by 
Eq. (2.2b) where the effective potential is evaluated at 
the minimum of the potential. 



P 



ai4ff , dV,s ) x' + PR 



-} 



(2.45) 



at 



dfj.R a^-" ^^^R J 2A 

The condensate density is defined to be pc = |(/)p 
the minimum of the effective potential. From (2.2c), the 
pressure is the negative of the effective potential 

2Xr 



P^-VR-x'\^f + 



1 



A* 



ix' + PR? 
4Afl, 

^ln[l 



r d3fc J- 1 

i (2^1 2 



Wfe - Efe - X + 



\Al 
(2.46) 



D. LOAF phase diagram 

As discussed in our previous papers, there are three 
possible solutions of the gap equations. These define 
three phase-space regions for a given coupling strength. 
They are as follows: 
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FIG. 1. (Color online) LOAF phase diagram. 



I. The "broken symmetry" case where </> ^ 0. From 
(2.44a), if (j) and <j)* are nonzero this means that 
x'^ — \A\'^ ■ Because of the U{1) invariance of 
the Lagrangian, we can choose and A to be 
real, in which case x' — A. In this region ojk — 
v/efc(efc + 2x')- Here < T < T,. 

II. The state when cj) = Q and < |A| < x' ■ In this 
region ujk - ^{tk + x'? -W-, for T, < T < T* . 

III. The normal state when (j) = 0, \A\ = 0, and x' > 0, 
for T* < T. In this region uJk = fife + x'- 

So the LOAF phase diagram shown in Fig. 1 is char- 
acterized by two special temperatures, Tc and T* . The 
critical temperature Tc corresponds to the appearance of 
the BEG condensate in the system. We will show next 
that the temperature T* corresponds to the onset of su- 
perfluidity in the system, and is related to the new order 
parameter, A. 



III. GOLDSTONE THEOREM 

The classical Lagrangian density with auxiliary fields 
and currents from Eq. (2.10) is given by 



L 



, \A{x)\'-x\x) 



2A 



J„(a;)$"(a;), (3.1) 



where is given in Eq. (2.12). With the exception of 
the currents, this Lagrangian is invariant under a global 
U{\) transformation of the form, 

0(x) ^ e*^ 0(x) , and A{x) e^'^ A{x) , (3.2) 

with x(x) unchanged. Consider the change in the La- 
grangian density under the infinitesimal change. 



$"(x) ^ $"(x)-H(5$"(x) , 



(3.3) 



From 



and using Lagrange's equation, we obtain 



- drp{x) + V • j{x) = - $„(x) g% JP{x) , (3.4) 



where 



p{x)^\4>{x)W 

i{x) = ^[cj,* {x)V^{x) - cj,{x) V0* (x) ' 



(3.5) 



Eq. (3.4) is a classical result and represents the 17(1) 
charge conservation equation. Here, we introduced the 
U{1) charge metric g°' p as the diagonal matrix given by 



g"^ =diag(l,-l,0,2,-2) 



(3.6) 



We multiply (3.4) by exp{ -S[<^, J] }, divide by Z, and 
integrate over the fields $ to derive a functional equation: 



i 1 

h'z\j 



yD$e-^[*"^l<D„(x)5"^ J^(a;) (3.7) 
^ { 1 d,p{x) + V • j(x)} . 



Integrating (3.7) over [da;] and discarding the surface 
terms gives 

j [dx]J^{x)g%^^ j D<i><i>„(x)e-^[*'^1 (3.8) 
= / [dx] J„(x)g";3$^[J](x) = 0. 



Changing functional variables from J to $ and using 
Eq. (2.35), we can write Eq. (3.8) as 

Differentiating (3.9) with respect to <i>"(x) gives 

= / |d.T'15;,'['f|(i,i')9",'f(a^') = -JuW/o ■ 
Expanding Q^^'" 0(x,x') in a Fourier series, 

g-^-''p{x,x') (3.11) 



P J (27r)3 



n— — oo 



and for constant fields and in the absence of sources, 
Eq. (3.10) gives [19], 



^-i"^,(0,0) #'3 = 0, 



(3.12) 
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where 



,0,2A,^2A* 
,0,-2A*,2A' 



(3.13) 



For the LOAF approximation, G^^'" 



-la 



/3(q,0) is given by 



-A* e, + x' 0* 
0* 7(9) S{q) S*{q) 

-cj, a{q) P{q) 

V -r Kl) 

(3.14) 



where 



a(g) =^-i'\[Xo](q,0) = ^-^'^ixoKq, 0) , (3.15) 

fi{q) = ^-i'i2[xo](q, 0) = ^-i'^[xo]*(q, 0) , 

7(g)-^-i'%[Xo](q,0), 

%) =^-i'°i[Xo](q,0) = ^-i'"2[xo]*(q,0) 

- ^"'''o[Xo]*(q, 0) - ^-i-'oIXoKq, 0) . (3.16) 

For small momentum values, the momentum expansion 
of the Green functions is 

a{q) = ao + aiq^ + • • • , (3-17) 
m = (io + Piq^ + • • • , 

7(9) = 70 + 71'?^ H , 

5{q) = 5o + 5iq^ + -- - . 

Since is missing the term proportional to x, we can 
extract the 4x4 matrix M"^, which is given by 



M"^(q,0) - A/o^^+Mf^g^^-.- 
Ze^+x' -A -r 

-A* eg + x' ~<l> 
a{q) /3(g) 
V -r f3*{q) a{q)y 



(3.18) 



where 

Mo"^(0,0) 

From Eq. (2.42), we find 
2n(/3wfc) + 1 



ao 



2uk 2€k 



(3.19) 



(3.20a) 



l3o 



(2^)3 AiJk { "^^^^^ + 1) 

+ 2n{Pojk) + 1 } , 

+ 2ri(/3wfc) + 1 1 . (3.20b) 



A^ f d^k 1 
~2 



Then the determinant of Mq is 

det[ Afo l^oP ) ( x" - I^P ) - 2ao x' |</>P 

+ /3o*A02+/3„A*0*2 + |(/.|4. (3.21) 

From the gap Eq. (2.44c) and for real fields we find that 

2A ■ 



(3.22) 



The Goldstone theorem corresponds to det[A/o] = 0. 
Hence, LOAF shows that Goldstone modes are present 
in region I and II of the phase diagram. The condition 
det[A/o] = is satisfied in region I because x' = 1^1, 
and 7^ 0, whereas in region II we have = and 
A ^ Q and the condition det[A/o] = is satisfied also. 
It can be shown that the determinant of the full 5x5 
inverse Green function also vanishes at q = and n = 0, 
det[^-i"/3(0,0)] = 0, when det[Afo] = 0. 



IV. SUPERFLUID DENSITY IN A BOSE GAS 
OF PARTICLES 

Traditionally, the density of the superfiuid component 
of the Bose gas is calculated using Landau's phenomeno- 
logical two-fiuid theory, as shown for instance in Ref. 20. 
In this theory, the superfiuid density is given by 



Ps = p - Pn 



(4.1) 



where p is the total density, ps is the superfiuid density, 
and the normal density, p„, defined as 



Pn 



67r2 



/•OO 

/ dkk* 


■ dn{pujk) ■ 


h 


dujk 



(4.2) 



where n{x) = [ 1 — e ^ ] ^ is the particle number density, 
and the dispersion relation is given by 



fife = 



2m 



(4.3) 



In this section, /3 = l/(fcBT). We will show next that the 
superfiuid density is proportional to the square of the 
order parameter A. To this end, we formally consider 
the series expansion of the normal density in powers of 
the order parameter, 

Pr^= p- kA^ + ■■■ , (4.4) 

SO that 

Ps = nA^ + --- (4.5) 
We begin by integrating Eq. (4.2) by parts. We use 
dn{j3u!k) i^k rn dn{j3ujk) m dn{/3ujk) 



dujk 



fife 



X fi^k dk 



ffik dk 



(4.6) 
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when A 0, to write Eq. (4.2) as 

_ _n^ [Zu. ^'^fc gn(/3^fc) 
G^^Jo ek + x' dk 

Integrating by parts gives 

1 /m\2 



(4.7) 



(4.8) 



dfc n 



{P^k) { 



4 41^1' 



')2 } ■ 



The second term in the last Hne of (4.8) is already of 
order , and in the limit A — >■ gives a term 



(2 



1 /m\2 r^^^nifiojok) f8 



i /my 



f-k + X ^3/ (efe + x') 



T^l^l • (4.9) 



For the first term, since Uk involves a term of order A^ , 
we need to expand it in a power series in A^: 



dn{PuJk) 



dA^ 



A' 



We use 

d{A^) = -d{ujl) = -2{ek + x'){^)dk. 
so that 

dn{f3ujk) 1 1 dn{(3uik) 



dA^ 2 VfiS'' + T^') dk 

this gives then for the first term, 



1 /TO- dn(j3ujk) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



2 V;i2^ fc(efe + x') afc 
p- k' A'^ + ■■■ , 



where k' is given by 



1 /my 



dfc- 



fc Wfc 



(2^)n;iy Jo {ek + x'r dk 

J_/TO\2 n(/3wofe) 



(4.14) 



(2t 



dfc- 
efe + X' 



1 



2£fc 

fife+X' 



where in the last line we have integrated by parts. Adding 
this result to what we found in Eq. (4.9) gives 



(27r)2 Jo Cfe + X' 

^ 2efc 8 el 
ek + X' 3 (cfc + X'? 



(4.15) 



This completes our derivation. 



V. MICROSCOPIC THEORY OF THE 
SUPERFLUID DENSITY 

Consider a superfluid moving with velocity v in the 
laboratory frame. The Lagrangian for this system is ob- 
tained by replacing the momentum operator for the sys- 
tem at rest by 



-Vi-^P=-V-?7iv. 

i i 



(5.1) 



Then, the superfluid mass density is given by the second- 
order derivative of the free energy. 



1 

V 



d^F{V,N,T, v) 

9t>2 



(5.2) 



where v is the velocity of the superfluid. The free energy 
is related to the grand potential ^1{V, /i, T, v) by 

F{V,N,T,v)^n{V,fi,T,v)+ pN . (5.3) 

It was shown in Ref. 17 that Eq. (5.2) is equivalent to 



1 

V 



dHiiv) 



11=0 



(5.4) 



which is what we use here. 

With respect to the laboratory frame, the Euclidian 
Lagrangian is given by 



C = 



1 



^ d(j){x) dcf)*{x) 



(5.5) 



where H is given by 



2m 



V V 



A 



(5.6) 



where the fields 4'{x) are given in Lab coordinates. In- 
troducing the Hubbard-Stratonovitch transformation by 
adding the auxiliary Lagrangian 



>Caux = -^ [x{x)-XV2\cf>{x)f^ 

+ l^\A{x)-X^'{x)\\ 



(5.7) 



to C gives 

C'=C + £aux 
_ 1 

~ 2 
where now 



^ d(f>ix) dcj)*{x) 

<P \x) (p{x) 



n' = (f>*ix) 



dr 



dr 



^ ^ 1 2 

— V • V H — mv 
i 2 



2to 

+ V2xix)-ti] Hx) 
-l[A{x) [rix)r + A*ix) [0(x)] 
+ ^[\A{x)f~x'ix)]. 



(5.8) 
(5.9) 

(5.10) 
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Adding currents, the action becomes 



1 



(5.11) 



dx] [dx'] M^) G-^\mx,x')cj)\x') 



+ J [dx] { [ \A{x)f - xHx)]/{2X) + .r{x) ^-{x) ] , 
where 

G-i^[$](x,x')=<5(x,x') {^^^^^^^^ , (5.12) 

with 

= K + d^, hi-^ = K-dr, (5.13) 
h^w^ h „ 1 2 fT. r . 

hy = rV-V + -mt; +V2x[x)-^. 

2m I 2 

Note that in the laboratory frame, the dependence of the 
action on the superfluid velocity is contained entirely in 
the inverse Green function G,7^"b[$] (x, x'). So comput- 
ing derivative of the action, we find 



dvi 



= - j[<^x\ji{x) , 



(5.14) 



where ji{x) is the classical mass current density. 



3^ {x)^-[<t>*{x)V,(t>{x)~(t){x)V, 4>* [x) ]-v,p{x), 

(5.15) 

with p{x) = |0(a;)p. The second derivative of S with 
respect to the superfluid velocity is simply 



dvi dvj 



Sij / [dx] p(x) 



(5.16) 



The partition function Z[J;V, p, P,v] is also a function 
of the velocity v. Then 



= / D<i> e~ 

zdv,. zj dv, 



[dx] {j,(x) 



(5.17) 



and 
1 d^Z 



Z dvi dvj 



0$ 



\ dvi ) V dvj J dvi dvj 



z 



j I [dx] [dx'] j D$ 
X - SijS{x,x') p{x) 

[dx] [dx']{j^(x)jj{x') - 6ijS{x,x') p{x)) , 



(5.18) 



However, we need derivatives with respect to the grand 
potential. These are given by 



dvi dv,. ' 



(5.19) 



d^Z 

dvi dvj 

So we find 



-/3Z 



dvi dvj 



-P^Z 



dfl dfl 

dvi dvj 



dvi 



1 dZ 

'JZ 'dv. 



[da;] (j',(a;); 



and 



1 d^Z 



dvi dvj dvi dvj pz dvi dvj 
1 dZ dZ 1 d^Z 



(5.20) 



(5.21) 



dv, dvj PZ dvi dvj 



[dx]{p{x)] 



[dx] [dx'][{j.{x) j,{x')) - {3,{x)) {3,[x))\ 



Assuming for simplicity that the superfluid is moving in 
the z direction, we obtain the superfluid density from 



Eq. (5.4), by evaluating Eq. (5.21) for i = j 
obtain 



Ps 



1 

V 



dHl{v) 
dv'? 



= P- Pn , 



where p is the total mass density given by 
1 



and pn is the normal mass density. 



1 

w 



[dx] [dx']jzzix,x') 



where now 

jij (x, x') = ( j, (x) jj (x') ) - ( j, (x) ) ( jj (x') ; 



We 



(5.22) 



(5.23) 



(5.24) 



(5.25) 



is to be evaluated at w = 0. That is, the action S now is 
evaluated at w = and the current is given by 

Mx) = ^[r{^)'^.<f>ix)-H^)^,rix)] , (5.26) 

without the additional Vi p{x) term. 

The calculation of the normal density from (5.24) re- 
quires the evaluation of the four-point correlation func- 
tion. 



(T.{r(xi)0''(x2)0^(x3)0'^(x4)}; 

1 5^Z[J] 



(5.27) 
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which in leading order in e is given by the products of 
two-point functions. 

G''\xi,X2)G'\x^,Xi) (5.28) 

+ G"^(xi, :r3) G^'(a;2, 2:4) + G''^{x2,x^) G'^'*(xi, X4) . 

In this approximation, the calculation of the current- 
current correlation function follows a straightforward 
path. We obtain 



Pn = jzz(0,0) 



(5.29) 



/•OO 


■ dn{l3ujk)- 


'0 





where jij (q, s) is the Fourier transform of the current- 
current correlation function in the sense of 

1 r +°° 
^-(^'^)-^/(^ ^ j.,(q,.)e^[-(--)-(--)l. 

^ s — — 00 

(5.30) 

Eq. (5.29) agrees with Landau's classical formula in 
Eq. (4.2). 



VI. JOSEPHSON RELATION 

The calculation of the superfluid density discussed in 
the previous section was performed in the laboratory 
frame. The dependence of the action on the superfluid 
velocity is contained entirely in the inverse Green func- 
tion, as discussed previously. This dependence can be 
removed by changing variables to a new field (/'t,(x), de- 
fined by 



(/.„(x) = e*(™^--)/''0(x) 



(6.1) 



which is sometimes called a "twist" transformation [21]. 
Then since 



-■i(mv-r) / h 



V — TO V 



(6.2) 



in terms of the fields 0i,(a;), the action (5.11) in the Lab 
frame becomes 



5[J,$„;At,/3,v; 
_ 1 
~ 2 



(6.3) 



[da;] [Ax'] 0„,(x) G-i%[$„](x, x') 0t(x') 
+ j [dx] { [ \A,{x)\^ - x\x) ]/(2A) + r{x) <J>^(x) } 
where now 



G-^\['^,]{x,x') = 5{x,x') 
with 



-Al{x) hi-) 



hi+) =h + dr, h^'^ =h-dr 



(6.4) 



(6.5) 



2to 



+ V2x{x) - A* : 



and we have defined the auxiliary field Ay [x) as 

Here the new fields and currents depend on the superfluid 
velocity v and are given by 

K^{^v,<i>:,X.A,,Al) , (6.7a) 
J«(a;) = (j„,j:,s,5„,5:), (6.7b) 

Jv a (2;) — ( Jt) I Jt) J , Sy ) , 



with similar definitions for the new currents Jy{x). The 
important thing to note here is that the entire depen- 
dence on V now resides in the fields rather that in the 
Green function operator. In particular, we have 

{^)^'i':{x)=mx,g"p<^^y{x), (6.8) 



where 



5% = (1,-1,0,2,-2) 



(6.9) 



The currents also now depend on v, but in the following 
derivation, we will set the currents to zero. So then 



(?) 



h\ dS 



h\ d^^{x) SS[<^>y 

[da;] TOa:^ .g-^ 



[a.,(«) 



dvi S^^{x 

SS[^y 



(6.10) 



and 



i I dvi dvj 



/ / [da;] [da;'] Xi a;' 



(6.11) 



-^'iix) Kix')\. 



d^^ix)S<i>^{x') 
So derivatives of the partition function are now given by 



\i) Zdv,^ Z J \i) dv,^ 



.12) 
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and 



h\-^ 1 d'^Z 



(?) 



i / Z dvi dvj 

( dvj ) dvi dvj 

rri^ J J [dx] [dx'] Xi x'j 



(6.13) 



(S$^^(r 



9% 9^. 



,5$^(x) 



(5$°(a;)5$g(x') 



-<i>(;(a:)<i>„^(x'))} 



Finally, using the results in Eq. (5.19), derivatives of the 
grand potential with respect to v is given by 



dVt _ I dZ 
dvi (iZ dvi 
im 



(6.14) 



hp J ^ ^ ^U$?(a;) > 



and 



/3 



dn dn 1 d'^z 



dvi dvj dvi dvj j3Z dvi dvj 
_ I dZ dZ 1 d^Z 
f3Z'^ dvi dvj PZ dvi dvj 



(6.15) 



+ 



JJ[dx] [dx'] x,x'^g"^gf\ 



5^^{x) S^^ix') 
The superfluid density is given by Eq. (5.4), 



Ps 



V 



■0=0 



(6.16) 



So we want to evaluate (6.15) at i ~ j ~ z and at v = 0. 
Now at zero currents. 



S'f'Six) 



(6.17) 



so the first term in the third line of (6.15) vanishes. So 
then the superfluid density is given by 



/ m 



\h^/3VJ 



) [[[dx][dx']zz' g^^^gP, (6.18) 



where the fields are now to be evaluated at v = 0. In 
leading order we can replace the fields by their 

expectation values so that the superfluid density in (6.18) 
becomes 



m 



[dx] [dx']z z'$„ ^-i^"^[$](x, x') $^ , 

(6.19) 



where are given in Eq. (3.13). Expanding 

Q~^''^ p{xtX') in a Fourier series according to Eq. (3.11), 
and inserting this into (6.19) and using the fact that the 
fields are constant, we obtain 



$a^"''";3[*](q,0)$'^ 



(6.20) 



Because the J7(l) charge of the auxiliary field x is zero, it 
is convenient to restrict ourselves to the (1, 2, 4, 5) set of 
indices. Therefore the relevant part of the inverse Green 
function G'^^'";3(0, 0) is the one we discussed previously 
in connection with the Goldstone theorem in Sec. Ill and 
is given by the matrix 



X 



-A* 








where we introduced the notations, 

a = ^-i'ii[xo](q,0) = ^-i'22[xo](q,0) 
/3 = ^-i-i2[xo](q,0) = ^-i'\[xo]*(q,0) . 



.21) 



(6.22) 



with D-^'^ j\x{^(x,x') given in Eq. (2.27). 

For small momentum values, the momentum expansion 
of the Green functions is 

a = ao + (^iq^ + • • • , (6.23) 
/? = /?o + + • • • . (6.24) 
Then, the inverse Green function can be written as 



e-i'"^[$](q,0) = AO + A./f^g2 
where 

-A 



P 



as before, and 



/ X' 

-A* 

\ 



X 





(6.25) 



(6.26) 



We find 



(h^jlm 

h^l2m 

ai /3i 

PI axi 



(6.27) 



-2x'|(/'P + 3(A0*2+A*(/)2) (6.28) 



8ao|A|2+4(A2/3*+A*2/3o) 
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and 

A/f ^ $^ -8 ai I + 4 (A^ ;3* + A* ^ . (6.29) 

We specialize now to the case of real fields and we recall 
that in LOAF we have (see Eq. 3.22) 



ao - /3o 



L(0,0)-i?-^-S(0,0) 



2A 



(6.30) 



For T <Tc, we have 7^ and x' = ^ > 0. Therefore, 
in region I, we obtain 

^2 8A^{ao-M, (6.31) 



which vanishes when we make use of Eq. (6.30). Simi- 
larly, ^aM^fj^^ vanishes for T > Tc, because (j) = 0. 
So, the Mqp contribution is identically zero everywhere. 
Therefore, we find that the superfluid density is given by 



Ps = 



n? 

fi2 



52 . 
lim— [I?- 

- jj- 



^(q,0)-^-i'i2(q,0) 
^(q,0) + i^-i'2^(q,0)] 



lim -^[^-i^ii(q,0)-^-i'i2(q,0)] . 

(6.32) 



In leading order, the superfluid density is 
f^m^ A^\ ^ , ^ . 



.33) 



In order to write (6.33) in terms of a component of the 
Green function, we compute the inverse of Eq. (6.21) and 
select the A4~^'*4 component. For real fields, this gives 

_ ia^-P')ix''-A^)-2ia-p)<j>^A 



1 



a(x'2 - A^) - X 

(6.34) 

For T <Tc where x' = ^ and when T > Tc where (/> = 0, 
(6.34) reduces to 2 (ai — /3i) in both cases. So we find 
that 

^ 2(ai-/3i), (6.35) 



lim 



q-^'o g2M-i'44(q) 
so that (6.33) can be written as 



Ps 



A ) 



lim 



1 



(6.36) 



A somewhat tedious but straightforward calculation 
shows that by inverting the full 5x5 inverse Green func- 
tion given in Eq. (3.14), we also obtain 

^—^2(a,-l3,)q^ + --- , (6.37) 

for both T < Tc and T > T^. So, Eq. (6.36) can be 
written as 

Ps — — ^ — lim — I (6.38) 

1~^U^Q^A{q.^) 

which is the Josephson relation. We note that the ex- 
istence of the pole in Q^Aiq^^), which is guaranteed by 
Goldstone's theorem when the U{1) symmetry is broken, 
is what leads to nonzero superfiuid density. 



VII. CONCLUSIONS 

To summarize, in this paper we derived the Joseph- 
son relation for a dilute Bose gas in the framework of an 
auxiliary-field resummation of the theory in terms of the 
normal- and anomalous-density condensates. In the self- 
consistent mean-field approximation to our auxiliary-filed 
formalism, the phase diagram of this theory features two 
critical temperatures, Tc < T*, associated with the pres- 
ence in the system of the BEG condensate and superfiuid 
state, respectively. For temperatures below Tc, a BEG 
condensate is present in the system, whereas T* signals 
the onset of superfiuidity in the system. As such, for all 
dilute systems of interacting Bose gases, the LOAF ap- 
proximation predicts the possibility of a superfiuid state 
in the absence of a BEG condensate in the temperature 
range between Tc and T* . The density of the superfiuid 
state is controlled by a second order parameter. The 
Josephson relation identifies this second order parameter 
as the square of the anomalous-density condensate. This 
result contradicts the usual Bose gas theory that does not 
feature an anomalous condensate, which predicts that the 
superfiuid density is proportional to the BEG condensate 
density. However, our theory is consistent with the case 
of dilute Fermi gases, where the BGS theory shows that 
the superfiuid density is proportional to the square of 
the gap parameter. In this sense, the auxiliary-field for- 
malism discussed here provides a unified approach to the 
study of fermonic and bosonic atom gases. 

Whereas the Goldstone theorem is a statement about 
the vanishing of the momentum-independent part of the 
determinant of the full inverse Green function of the the- 
ory, the Josephson relation is related to the coefficient of 
the determinant proportional to the square of the mo- 
mentum. The Josephson relation for the superfiuid den- 
sity reads 



8to2 A^ 



lim 



(7.1) 



where G'*yi(q, 0) is the propagator corresponding to the 
expectation value {A* A), evaluated for zero energy 
transfer. These results will remain preserved at all or- 
ders in the many-body approximations of the auxiliary- 
field formalism. 
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